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Abstract 

A source term in the quantum Boltzmann equation, which accounts for the 
spontaneous creation of e"'"e~-pairs in external electric fields, is derived from 
first principles and evaluated numerically. Careful analysis of time scales 
reveals that this source term is generally non-Markovian. This implies in 
particular that there may be temporary violations of the if-theorem. 
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The evolution of the quark-gluon plasma, believed to be formed in the 
course of relativistic heavy-ion collisions, is commonly described by means 
of a transport equation |1|, |^, |, |]. It is well understood how a trans- 
port equation can account for acceleration in external fields, scattering, or 
(hadro-) chemical reactions of the microscopic constituents. There is, how- 
ever, another physical process which becomes increasingly important at high 
energies: regions of very large chromoelectric field strength may develop and 
subsequently decay by emitting quark- ant iquark pairs ^ . This gives rise 
to the fragmentation of chromoelectric flux tubes ('strings'), a mechanism 
frequently invoked to model hadron production [^, ^, How such spon- 
taneous creation of particles can be incorporated into a transport equation 
is still not fully understood. 

Clearly, the transport equation has to be modified by a source term. What 
is this source term? How can it be derived from the underlying microscopic 
dynamics? These issues have recently been approached in a Wigner func- 



tion formulation [11, |T^, 13, 0, 15 1. But, aside from the fact that it lacks 
an intuitive probabilistic interpretation, this approach suffers from several 
practical limitations. The source term cannot be determined completely: it 
is not known how the longitudinal momenta of the produced particles are 
distributed. It has been suggested that the distribution is a J-function M; 



but while such an ansatz may be useful for practical purposes [0, |T^, |T5 
it is certainly not exact. Furthermore, an interplay of pair creation and col- 
lisions — possibly leading to a modification of the source term — has not 
yet been considered. And finally, the Wigner description is not suited for 
discussing the apparent irreversibility of the particle creation process or the 
associated generation of entropy. Since pair creation in an external field is 
merely a single-particle problem (see below), the Wigner function retains 
complete information about the microscopic state of the system. Yet irre- 
versibility never manifests itself on the microscopic level; it only emerges 
after a suitable coarse-graining. 

I choose a different approach. In collision experiments one usually mea- 
sures the momentum distribution of the outgoing particles; i.e., one deter- 
mines the occupation n^^^p, t) := {N^{'p)){t) of the various momentum states, 
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with the number operators given by 



N+{p) ■= '^b\p,m^)b{p,m^) . (1) 

(Here a and b denote particle and antiparticle field operators, respectively, 
p the momentum and the spin component.) This suggests attempting 
to describe the evolution of the occupation numbers {n:f{p,t)} directly and 
to derive a kinetic equation for them — including the source term — from 
first principles. I will do so with the help of a very powerful and broadly 
applicable tool: the so-called projection method. This method, pioneered by 



Nakajima Zwanzig |T8|, |T9[ and others 0, |21|, |2|, |2|, ||, |5[, is based 
on projecting the motion of the quantum system onto a low-dimensional 
subspace (the 'level of description') of the space of observables (Liouville 
space). It allows for a clear definition of crucial concepts like the memory time 
or the coarse-grained entropy, making it especially suited for an investigation 
of the irreversible features of the dynamics. 

As in the works cited above, my investigation is based on a simple model 
from quantum electrodynamics. I consider the creation of e'^e'-pairs in a 
homogeneous, time-independent electric field E, a process often referred to 
as the Schwinger mechanism The starting point is the Dirac 

equation 

zn^^mt)) = H\m) (2) 

with 

H = ^p-d + ml3 + qAo (3) 
and Ao(r) = —E ■ f {q = —\e\ for electrons). We further define p{t) : = 



p + qEt, the transverse energy e± := Jm'^+pj_, the total kinetic energy 



e[p(t)] := ^^Je']_ +p\\{ty, and the dynamical phase 

<Pf^■.= -l£'dt'e\p{t')] . (4) 



'Longitudinal' and 'transverse' refer to the direction of the electric field. 
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The eigenstates |i, ±) = |p(tj), m^, ±) which correspond to momentum 
p{ti), spin component and positive or negative energy ±e[p(tj)], evolve 
according to 





The evolution thus mixes positive and negative energy eigenstates, with 
respective amplitudes afi and Pfi] equals the probability for having 

created an e"'"e~-pair with (final) momenta ±p{tf) during the time interval 
[ti,tf]. The amplitudes are determined by the differential equation 

(6) 

with initial conditions an = 1 and Pa = 0, and the dot indicating differenti- 
ation with respect to t/. 

In view of applying the projection method, the above results have to be 
translated into the language of field operators. To do so, I will use the for- 
mulation of quantum statistical mechanics in Liouville space [^. There 
the evolution of (Heisenberg picture) operators is determined by the so- 
called super-operators £ ('Liouvillian') and U; these super-operators play 
a role analogous to that of H and U in Hilbert space. Employing the short- 
hand notations aj = a{p{tj),mz) and b-j = b{—p{tj),—mz) for the par- 
ticle and antiparticle field operators, and making use of the general rule 
U{t2,ti)a^ ijp) = {U(t2,ti)ip), one finds 

the evolution law for (a, b^) follows by Hermitian conjugation. Thus pair 
creation can be described by a time- dependent Bogoliubov transformation 



The Liouvillian 

dt2 



UihM) (8) 

«2=tl 



may be written as the sum 

C = £diag + SC (9) 
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of a diagonal part, responsible for acceleration, and an off-diagonal part which 
is responsible for the mixture of particle and antiparticle states, i.e., for pair 
creation. With the definition Pu := (321 , the latter is given by 

t2=tl 

Starting from the above microscopic equations, we now want to derive a 
kinetic equation for the occupation numbers n:i^{'p,t). Provided the initial 
state is the vacuum, 

P(io) = |0)(0| , (11) 

momentum and charge conservation dictate n^{p,t) = n^{—p,t) for all later 
times t; it then suffices to consider the evolution of only, say, the electron 
occupation numbers n-{p,t). Their evolution equation must have the struc- 
ture 

n_ {p, t) + qE- Vp-n_ {p, t) = ri!_°" {p, t) (12) 

with some source term rf°'^{p,t). Since we know that this source term ac- 
counts for transitions between positive and negative energy eigenstates, it is 
tempting to write down a rate equation of the form 

h'°^{p, +e, t) = ir(p ) ■ [n{p, -e, t) - n{p, +e, t)] , (13) 

r{p) being the respective transition rate. If this were correct, the identifica- 
tions n{p,+e,t) = n^{p,t) and n{p,—e,t) = 2 — n^{—p,t) would then lead 
to 

h^--{p,t) = r{p)-S{p,t) (14) 

with 

S{p,t):=l-ln4p,t)-ln4-p,t) . (15) 

Such a source term, however, can only be correct in the Markovian limit 
— an approximation which is not always justified. Careful investigation 



31, 321 reveals that the above ansatz for the source term has to be modified: 



assuming the quasistationary limit (to —oo) one finds 

roc 

\m= dTR{p,T)-S{p-qET,t-T) . (16) 
Jo 
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This source term involves an integration over the entire history of the system, 
thus accounting for finite memory effects and rendering the evolution of the 
occupation numbers generally non-Markovian. 

The kernel R{p, r) can be obtained with the help of the projection method 



[3^ . One key ingredient in the derivation is the introduction of a super- 



operator Q which projects onto the irrelevant degrees of freedom; in our 
case, 

QN^ = QN+ = , (17) 
whereas other combinations of field operators are unaffected: 

With this definition one finds 

R{p,t) = -{0\SCexp{iQCQT)S£N_{p)\0) , (19) 

a general result which holds for arbitrary field strengths. 

The formal expression for the source term can be easily evaluated in the 
limit of weak fields. Provided E -C m'^/hq, then |/?| <^ |0| and hence 6C may 
be regarded as a small perturbation. In this case it is legitimate to replace 

expizQCQr) exp(?Q/:diagQr) , (20) 

leading to 

/•CO . ... 

n^_°"(p,t) = 4Re / dr P*{-t, -r) e''^'^^-^'^^ ^(0,0) S{p- qET,t - r) (21) 
Jo 

(with j3{t2,ti) = P21 and 0(t2,^i) = 02i)- This source term is consistent 
with the Schwinger formula: assuming that the system is dilute {S = 1), 
using the differential equation (H) in the limit ~ 1, and employing the 
Landau- Zener formula ||3^, |3^, one finds that 



w :- 



a result which does indeed agree with the leading term in the Schwinger 
formula. 
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The source term (|2T| ) may be evaluated numerically. For simplicity I will 
take the system to be dilute, 5" = 1; the source term then no longer depends 
on t. It is convenient to introduce 

a := hqE/el (23) 

and to consider, rather than the source term itself, the dimensionless quantity 




The pre-factor in its definition has been chosen such that for = 0, ?7 is of 
the order one. As the pre-factor is independent of p\\, t] will correctly describe 
the distribution of the longitudinal momenta of the produced electrons. One 
can show that 

^^"'^"/^-^ = (^) 2[1 + ip,/e,r] /-l^"coshV(L||A.) (^) 

(25) 

with ip defined implicitly as the solution of the equation 

sinh (/9 cosh (/} + Lp = x + {p\\/e±)^Jl + {p\\/t^_Y + sinh~-^(p||/e^) . (26) 

I calculated t] numerically, using a combination of Filon's integration for- 
mula with an efficient root-finding algorithm. The results for weak fields 
(a < 1; figure 1) may at first seem surprising. Clearly the momentum dis- 
tribution of the produced electrons is not narrowly peaked around p\\ = 0; 
it is neither a (5-function nor a thermal distribution. Rather, electrons are 
being produced predominantly in the direction of the external field {p\\ > 0). 
Electrons moving in the opposite direction (py < 0) are being annihilated: 
for them, the production rate is negative. Of course, such negative produc- 
tion rates are sensible only if there are electrons available for annihilation. In 
the quasistationary limit this is the case: electrons which have been emitted 
with positive momentum are subsequently being decelerated and may then, 
as soon as p\\ < 0, be (partly) annihilated again; there remains a small sur- 
plus which manifests itself as a positive total production rate. As another 
surprising feature, t] displays (approximately) periodic oscillations whose pe- 
riod scales with a. This may be understood qualitatively if one views pair 
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creation as a tunnelling process from the negative to the positive energy 
continuum The barrier between these continua has a spatial width of 
the order e^/qE, inducing a 'momentum quantisation' Apy ~ hqE/e± and 
thus A(p||/ex) ~ a. Interference of multiply reflected electron wave func- 
tions then leads to the observed oscillations. For a strong field (a > 1; figure 
2), the naive tunnelling picture breaks down; both the oscillations and the 
annihilation of particles (negative rates) become less pronounced. 

As we discussed previously, the source term is generally non-Markovian. 
It exhibits two characteristic time scales: (i) the memory time rmem(p), which 
corresponds to the temporal extent of each individual creation process and 
which indicates how far back into the past one has to reach in order to predict 
future occupation numbers; and (ii) the production interval rprod(p) — the 
inverse of the production rate — , which corresponds to the average time 
that elapses between creation processes and thus constitutes the typical time 
scale on which the occupation numbers change. Only if r^em ^ Tprod can 
memory effects be neglected and the evolution be considered approximately 
Markovian. 

In the weak field limit both time scales can be extracted from the source 
term (pT]). First the memory time: The factor 

7 / r,,2 (27) 

{T-pii/qEy + {e±/qEy 

constitutes a Lorentz distribution in r, centered around p\\/qE with width 
e±/qE. Significant contributions to the source term thus come from times r 
which are smaller than (py + e±)/qE. As the typical momentum scale is set 
by Aj9|| ~ hqE/e±, we may conclude 

Tmem ~ \ ^ ■ (28) 

e_L qE 

The memory time combines two time scales of different origin, (i) The time 
h/e± is proportional to h and therefore of quantum mechanical origin. It 
corresponds (via the time-energy uncertainty relation) to the time needed to 
create a virtual particle-antiparticle pair, and may thus be regarded as the 
'time between two production attempts.' (ii) The time e±/qE, on the other 
hand, is independent of h and therefore classical. It can be interpreted in 
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various ways, depending on the picture employed to visualize the pair cre- 
ation process. If pair creation is viewed as a tunnelling process, the classical 
memory time coincides with the time needed for the wave function to tra- 
verse the barrier with the speed of light Alternatively, pair creation may 
be viewed as a non-adiabatic transition between the two time-dependent en- 
ergy levels ±e[p(t)]. In that case the classical memory time corresponds to 
the width of the transition region, i.e., the region of closest approach of the 
two levels. Finding the production interval is less straightforward. Assuming 
p\\ = for simplicity, again invoking the weak field limit, and exploiting the 
fact that the source term must be consistent with the Schwinger formula, 
one can show that 

V„.(0.pj~i|exp(^) . (29) 



As long as E <^ m? /hq < e\/hq, the particle creation process is Marko- 
vian: Tmem ^ Tprod- In the weak field limit, therefore, the entropy associated 
with the coarse-grained level of description spanned by {N±{j))}, 



n_{p,t) ^^ n_{p,t) ^ / ^'^^'^"^ ^Ind ^^^P^^'> \ + 



+ (ri„ ^ n+) , (30) 

obeys an if-theorem. {Q denotes the volume.) The monotonous increase 
of the coarse-grained entropy explains why spontaneous pair creation is per- 
ceived as irreversible. This apparent irreversibility is, of course, a conse- 
quence of the coarse-graining: information is being transferred from acces- 
sible (slow) to inaccessible (fast) degrees of freedom. The slow degrees of 
freedom are the occupation numbers of the various momentum states. From 
these, information gradually 'leaks' into unobserved degrees of freedom: cor- 
relations and rapidly oscillating phases which entangle the respective wave 
functions of the members of a particle-antiparticle pair. 

As soon as -E > m'^/hq, the situation changes. The source term (pT]), which 
was derived in the weak field limit, is then only a rough estimate. Already this 
weak-field estimate becomes non-Markovian: at E = rn^/hq the production 
interval and the memory time are of the same order r ~ h/m. This is a 
clear indication that at this point conventional Markovian transport theories 
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must break down. There may be temporary violations of the if-theorem: the 
coarse-grained entropy, while still increasing on average, may now oscillate 
(on the same scale r ~ h/m). Such oscillations have indeed been observed 
in numerical simulations []T3[. A systematic study of these memory effects 
should proceed from the general equations (|16D and ([T9|) . Although such 
an enterprise is beyond the scope of this letter, we can already say that 
(i) memory effects become significant at large field strengths; and (ii) the 
projection method can account for these memory effects and thus appears to 
be a suitable tool for their investigation. 

The above analysis can be extended to include binary collisions of the 
produced particles. This is done by replacing 6C 6C + V, where V contains 
the two-body interaction. To lowest order perturbation theory, pair creation 



and collisions do not interfere ||31|; the additional interaction gives rise to a 
separate collision term. Like the source term, this collision term is generally 
non-Markovian and must be subjected to a time scale analysis, leading again 
to a criterion for the validity of the Markovian approximation. One finds that 
there are two contributions to the memory time: the average time needed 
for a particle to pass through an interaction range, and the typical 'off- 
shell' time given by the time-energy uncertainty relation. For the Markovian 
approximation to be valid, these have to be smaller than the average time 
that elapses between two successive collisions pi], p7 . 



Let me summarize the main conclusions, (i) The source term in the quan- 
tum Boltzmann equation can be derived in an unambiguous fashion by em- 
ploying the projection method, (ii) To lowest order, the source term is not 
altered by the presence of collisions, (iii) In the weak field limit, E -C m? /hq, 
the source term is given by (^) or (pSf) , respectively. It is then Markovian, 
and the coarse-grained entropy increases monotonically. As information is 
continuously being transferred to inaccessible degrees of freedom, sponta- 
neous pair creation appears irreversible, (iv) But as soon as > m?/hq, 
there may be sizeable memory effects, leading to temporary violations of the 
if-theorem. Their description is beyond the scope of conventional Marko- 
vian transport theories. A more suitable starting point appears to be the 
projection method, in particular equations ([T6| ) and (p^Q]). 
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Figure Captions 



Figure 1: The re-scaled production rate 77 as a function of p\\/e±_ for weak 
fields (a = 0.2, 0.3, 0.7). 

Figure 2: The re-scaled production rate 77 as a function oip\\/e±_ for a strong 
field (a = 2.9). 
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